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A theoretical model is developed describing high-order harmonic generation (HHG) from a gas 
of multiply-charged ions driven by a laser field of relativistic intensity. Macroscopic propagation of 
harmonics is investigated in a relativistic HHG setup where the relativistic drift is suppressed by 
means of x-ray field assistance of the driving laser field. The possibility of phase-matched emission 
of the harmonics is shown. The laser field geometry is optimized to maximize the HHG yield with 
the corresponding phase-matching schemes. Crucial issues determining the macroscopic HHG yield 
are discussed in detail. 
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I. INTRODUCTION 

High-order harmonic generation (HHG) is a reliable 
source of coherent soft x-ray radiation in the nonrela- 
tivistic regime. With the state-of-the-art technology, co- 
herent x-rays of up to ~ 3 keV pQ of photon energy 
can be generated. The most favorable conversion effi- 
ciency for nonrelativistic multi-keV harmonics is antici- 
pated with mid-infrared driving laser fields at high gas 
pressures j2rH] ■ A further increase of the photon energy 
can in principle be achieved by increasing the laser in- 
tensity. However, the applicable laser intensity is limited 
in two ways [5j. First, the relativistic electron drift pre- 
vents recollision and results in a dramatic suppression of 
the HHG efficiency. And second, the strong field causes 
rapid ionization of the medium leading to a large free 
electron dispersion and along with a significant phase- 
mismatch. The electron recollision is suppressed when 
the drift distance becomes larger than the electron wave 
packet size at the moment of recollision [5] . This happens 
when the laser intensity exceeds 10 16 — 10 17 W/cm 2 for 
infrared (IR) wavelengths. The ponderomotive potential 
of the laser field in this case amounts to U p ~ 3 keV and 
the achievable cutoff frequency for HHG to lo c 10 keV. 
This indicates the limit of nonrelativistic HHG. 

Various methods to counteract the relativistic drift 
have been proposed. To suppress the drift, different laser 
field geometries, in some cases with an additional field, 
can be applied [7H20). Highly charged ions moving rela- 
tivistically [2TJ 122] or a gas of positronium atoms [231 HI] 
can also be employed for this purpose. However, all these 
efforts have only addressed the drift suppression problem 
for the emission from a single atom rather than coher- 
ent emission from a macroscopic gas target where phase- 
matching becomes crucial. For the first time both prob- 
lems of relativistic HHG, namely, the relativistic drift 
and the phase-matching, have been solved at the same 
time in [35] where the macroscopic yield of HHG has 
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been calculated in the setup consisting of two counter- 
propagating attosecond pulse trains [TTJ]. In this setup, 
the relativistic drift caused by the ionizing laser pulse is 
reverted by the counter-propagating pulse inducing re- 
combination. It appears that, specific to this setup, an 
additional harmonic phase exists which depends on the 
time delay between the driving pulse trains. We have 
shown in [25] that this additional phase of the emitted 
harmonics can be tuned to compensate the phase mis- 
match caused by the free electron background. How- 
ever, the setup of counter-propagating pulses is rather 
challenging, for instance, the requirement for a small 
pulse distortion imposes a rather strong restriction on 
the medium length. Additionally, a precise modulation 
of the laser intensity along the propagation direction is 
required. 

Another appealing scheme for relativistic HHG ex- 
ists based on XUV assistance [TH] which seems experi- 
mentally less demanding than the scheme with counter- 
propagating attosecond pulses. The usefulness of XUV 
light assisting a strong laser field has been demonstrated 
in the nonrelativistic regime for various purposes. It has 
been used to enhance HHG by many orders of magnitude 
compared with the case via a fundamental laser pulse 
alone [13 [27j. When the XUV field has the form of an 
attosecond pulse train a single quantum path can be se- 
lected to contribute to HHG and in this way allowing 
to manipulate the time-frequency properties of harmon- 
ics as well as to enhance a selected bandwidth of har- 
monics [28 30J. Tuning the XUV field to a resonance 
between a core and valence state can lead to the emer- 
gence of a second plateau that is shifted to higher energies 
by the former resonance energy with respect to the first 
plateau [31]. 

In the relativistic regime the XUV/x-ray assistance 
can be employed to overcome the relativistic drift mo- 
tion [15] . Thereby, the XUV frequency requires to ex- 
ceed the ionization energy to liberate the electron with 
a single photon and to deliver a significant initial mo- 
mentum to the freed electron. This way, the electron 
can obtain sufficient momentum in the direction oppo- 
site to the laser propagation direction to compensate for 
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subsequent drift motion and return to the atomic core, 
recombine and emit harmonics after the excursion in the 
relativistically strong laser field. The medium is a gas 
of multiply charged ions with an ionization energy large 
enough to withstand the strong optical laser field. How 
much the XUV assisted setup for relativistic HHG fa- 
vors phase- matching needs investigation. Formalisms de- 
scribing macroscopic effects due to ionization and phase 
matching in HHG are restricted to the non-relativistic 
regime, for a recent review see, e.g., [32] . 

In this paper, we investigate the feasibility of phase- 
matched emission and the macroscopic yield of harmon- 
ics in the relativistic regime of the x-ray assisted HHG 
setup in a strong IR laser field. Generally, the efficiency 
of HHG is rather small even in the nonrelativistic regime 
due to the wave packet spreading. In the relativistic 
regime, the single-atom HHG emission rate continues 
to decrease even when the relativistic drift is compen- 
sated [25| . Thus, a large phase-matching volume is cru- 
cial in order to achieve a significant HHG yield. Fur- 
thermore, the large ponderomotive potential is likely to 
result in rapid phase changes if ions emit under different 
conditions. For generating relativistic harmonics both 
challenges have to be met: circumventing the drift and 
having the setup stable against phase changes. Our pre- 
sented setup overcomes both issues and renders a mea- 
surable HHG yield in the relativistic regime possible. 

The structure of the paper is the following. In Sec. [Hj 
the theory of macroscopic HHG is presented applica- 
ble for any field geometry in the relativistic regime. In 



Sec. Ill the developed theory is applied to calculate the 



macroscopic HHG yield for the setup of x-ray assisted 
relativistic HHG. Our conclusion is presented in Sec. [TV] 



II. MACROSCOPIC MODEL FOR 
RELATIVISTIC HHG 

A. Macroscopic HHG yield 

In this section, our model is presented for the cal- 
culation of the harmonic spectrum from a macroscopic 
gas target suitable for relativistic laser intensities. In 
the non-relativistic regime, the standard approach for 
the calculation of the macroscopic HHG response in- 
corporates the single-atom contribution via the time- 
dependent dipole moment [32H3I]- However, employing 
the dipole moment for the radiation response assumes 
that an emitted harmonic wavelength is much longer than 
the spatial extensions of the emitter. This approach fails 
for sufficiently small wavelengths because then the re- 
tardation between different points of the emitting wave 
packet becomes important [35] [35] . Our approach uses 
the complete current density distribution of each atom 
rather than the dipole moment. Retardation between 
different emission points within the distribution is taken 
into account by a phase factor. The link between the mi- 
croscopic (atomic) current density j and the macroscopi- 



cally emitted harmonic electric field Eh is obtained from 
Maxwell's equations similar to the non-relativistic ap- 
proaches [3"2"HMj . The Fourier component of the emitted 
harmonic electric field from a gas target is given by [37] 



E h (x',<j h ) = i 



R 



where j(x, wh) is the Fourier component of the current 
density, wh is the frequency of the emitted harmonic 
light, k.H = fcnn' the wave vector, R = |x — x'| the 
distance between the emission and observation point, x' 
the coordinate of the observation point and n' = x'/lx'| 
the unit vector in the observation direction (see Fig. lj. 
Absorption of the harmonic photons is neglected [53] be- 
cause their energy is much higher than the largest atomic 
transition energy. The current density j is exclusively 
determined by the HHG process. For the evaluation of 
Eq. ([!]), we restrict ourselves to the far-field zone which 
is sufficient for calculating the overall HHG photon yield. 
The far-field zone is determined by the conditions that 
the distances from the emitters to the observation point 
are larger than the wavelength of the emitted radiation 
as well as the size a of the emitting region (k^R 3> 1 
and R ^> a). We thus can expand Eq. (fil) over a small 
parameter a/R using R = |x — x'| ~ |xf — x • n' |37j . 
When inserting this expression into Eq. ([I]), the expo- 
nential function splits up into two parts. One term con- 
tains x' and is thus a general phase factor depending 
on the constant observation point that can be separated: 
E h (x',wh) = e lfeH l x 'lE H ,o(n', wh). In the following we 
consider Eh.o(h',wh) only and find 

E H ,o(n> H ) =i ||/ d3xn ' x (J(x,«h) x n>- ikH X . 

(2) 

The total current density distribution consists of a sum of 
the current densities of the single atoms j a (x a , x, t) with 
positions x a : 



j(x,w H ) 



d 3 x a p(x a ) / di j a (x a , x, t)e 



(3) 



where p(x a ) is the atomic number density. Inserting 



laser pulse 




yield 



FIG. 1: (color online). Geometry of the medium and detector 
including the definitions of the coordinates. The dashed lines 
denote the divergence angle of the harmonic radiation. The 
box on the right schematically shows the measured angular 
distribution. 
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Eq. ([3]) into Eq. ^ yields the final expression for the 
macroscopically emitted harmonic field 

E H ,o (n',UH) 

= / d 3 X a ,(x a ) / d 3 X 

y"dtn' x (j a (x a ,x,t) X n ') e -^HX.n'+ 1WH * 

= J d 3 x a p(x a ) n' x (j a (x a , n', w H ) X n') (4) 

where j a (x a , n', w H ) = / d 3 x / dij a (x a , x, t) exp{-ifc H x • 
n' + iuj^t}. Note that the combination of outer products 
in Eq. Q prevents mathematically from emission in the 
direction of the current density vector. However, due 
to phase-matching, the macroscopic emission is in many 
cases mainly along the propagation direction of the laser 
n' « z, (see Fig. [I]). When the laser is linearly polarized 
in x direction, the current density vector is parallel to the 
x — z plane. Thus, we can approximate n' x (j x n') = 
j — (n' • j)n' rj j x x and can restrict ourselves to the x 
component 

^H,o,x(n', w H ) = i-^| J d 3 x a p(x a )j aiX (x a ,n',a;H) (5) 

to describe emission in this case. 

The overall emitted energy can be obtained via inte- 
grating the Poynting vector S(r, t) = ~Eg(x', t)n' over- 
all emission directions in the far field 



W 



cR 2 

4<7T 



dt J dfiE^(x,t) 



(6) 



By inserting 



E H (x',i) = — / dw H e- iWHt EH(n',u;H) (7) 

into Eq. the energy can be calculated via an inte- 
gration over the spectrum [55] : 



W 



cR 2 
(2tt) 2 



dn / dw H |EH,o(n',w H ) 



(8) 



The emitted spectral photon number per solid angle from 
Eq. (|| is given by: 



dN 



cR 2 



dcjjjdri (27t) 2 u;h 



|E H ,o(n',WH)| : 



(9) 



B. Single-atom current density 



in the relativistic regime via the Klein-Gordon equation. 
The Klein-Gordon current density of a particular atom 
at position x a in the laser field is given by 

j a (x a ,x,t) = (*;>)j*x>) + **.(a:)j**£>)) (10) 

where j = p + Al(x)/c and ^ Xa (x) is the solution of 
the Klein-Gordon equation when the binding potential 
is centered around x a . In the following, the time-space 
coordinate is x = (ct,x), the wave four- vector of the 
laser field = (wl/c, Isl), and the metric tensor = 
diag(l, -1,-1,-1). 



By a Fourier transformation of Eq. ( 10 ) and partial 



integration, the spectral current is obtained 
j a (x a ,cj H ,n') 



dW fc H*(v^>)j* XE >) + * x >)j***>)) 



d*xe ik * x V* (x)(j-ik H )*x a 0r) 



(11) 



We calculate the electron wave function in the field of 
the laser and the ionic core by means of the strong-field 
approximation (SFA) [101 HI] 

9^(x) = <Hx-x a ,<)+ / d 4 x'G Ij (x,x')Vi(x')ct>(x'-x a „t') 

(12) 

where V\ — 2i(A L (7j)/c)V — A L (?7)/c 2 ) is the term in 
the Hamiltonian describing the electron interaction with 
the laser field with the vector-potential A l (77) and phase 
t] = kh ■ x. The Volkov propagator Gt,(x,x') in a plane 
wave laser field is given by [121 E3] 



G h (x,x') = -i0(t-t') 



id 3 q 



exp [-iSi,(x,x')\ 
(13) 

with the classical action of an electron in the laser field 



2e q (2vr) ; 



S hM = - ,0+ /\ [ (q + A - ( f 2c) - AL( ^ )/c 

J,f L fc L • q 

(14) 

the energy-momentum four- vector q = (e q /c, q), 
the energy e q = ^/c 2 q 2 + c 4 . Here <f>(x) — 

j^i-h) exp{ ~ i[(c2 - j p }t + x • Al/c]} ' where ^° (x) 

is the nonrelativistic ground-state wave function. Insert- 
ing Eq. \12\ into Eq. ( 11 ) and applying the usual assump- 
tions |43U44j of neglecting bound-bound and continuum- 
continuum transitions and the time-inverted process, we 
obtain 



In Sec. |II A[ the macroscopic HHG yield has been cal- 
culated via classical electrodynamics. Since the expres- 
sion Eq. |5| for the emission field relies on the current 
densities of a single atom in the gas, we continue to derive 
the single-atom current density quantum mechanically 



j a (x a ,cj H ,n') 

2 f d 4 x I dVe ifcH -^*(x-x a ,t)(j- lk H ) 



(15) 
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Electron wave function in a distorted plane 
wave laser field 



The Volkov propagator Eq. |13j) describes the evolu- 
tion of the wave function of an electron in a plane wave 
laser field with vector potential Al(?7). However, in many 
practical situations this assumption on the laser field is 
not met, in particular, when a focused laser field or mul- 
tiple laser beams are applied [4"5H4"7] , or when the disper- 
sion distorts the laser pulse. For this reason, we need to 
find the electron wave function in an external laser field 
where the vector potential depends not only on phase 77 
but also on the position z along propagation. The devia- 
tion of the laser field from the plane wave form is assumed 
to be a perturbation, so that the total vector potential 
reads 



A(t7,x) = A L (7y) + A P (?7,x), 



(16) 



with |Ap(77, x)| <C |Al(??)|. We find the solution of the 



Klein-Gordon equation in the field of Eq. ( 16 ) using the 



cikonal approximation, see e.g. j48J[49], in which the im- 
pact of the perturbation onto the wave function is taken 
into account by an expansion of the wave function phase. 
The electron wave function is described by the Klein- 
Gordon equation 



(<9^ + c 2 )^(x) = V^(x) , 



(17) 



with V = 2i(A(?7)/c)V - A 2 (r))/c 2 . In order to solve 
Eq. (17), the ansatz 



= / nln ^ exp{-i[S L (x, x Q ) + S P (x, x )]} . 



(18) 



is employed. Here, St,(x, xq) is defined in Eq. (14 1, 
with an arbitrary constant xq, and, consequently, 
exp(— iSi,(x, xq)) solves the unperturbed equation (17) 
with A(j], z) = Al(?7). Inserting the ansatz Eq. (18) into 
Eq. fl7|), one finds 



1 



-id 2 S P - 2d t S L d t S P - (d t Sp) 2 } 
-[-iV 2 S P - 2VS , l VS'p - (VSp) 2 ] 



c 

A 2 



+ 2 — V(S L + S P ) 
c 

A L A P 



(19) 



In the applied approximation, |V5 P | -C |Ap/c| and 
\d t Sp\ < |Ap|, which allows to neglect the (VSp) 2 
and (d t Sp) 2 terms. When additionally £ = El/cujl ^> 
uj l /c 2 , one has |V 2 S p | < \(A P /c)VS p \ and \d 2 S p \ < 
\dtSt,dtS p \. Then, V 2 <Sp, d 2 Sp terms also can be ne- 
glected, yielding 

n A L Ap 



4 2 
Ap 



d t S L d t Sp + 



2^S L 

c 

-2VS L + — (Al 

c 



(20) 



The equations for the characteristics of the first order 



partial differential equation ( 20 1 are 
dt{u) 



du 
dr(u) 
du 



= -2V 1 5 L + -(A L +A P ). 
c 



It follows from Eq. (|21) and Eq. (|22j that 

dr dt 

- 2fc L -p = k L - ^l^- 

ou du 



and 



-(k L r - iu L t)/(2k L ■ p) = r,/(2k L ■ p) 



{21) 
(22) 

(23) 
(24) 



Integrating Eq. ( 20 ) and employing Eq. ( 24 ) we derive 



[ v 1 

Sp{x,x )= / dry- 



VSt 



A h (fj) Ap(?j,x L (x,^rj)) 



2c 



xA P (?7,x L (x,?7,77))/c 



(25) 



with 



x l(x, 77,77) = x- / dr]'—^— \] 
Jr, kh-p[ 



Al(t?') 



kh -P 



A h (r,'^A L (ri') 
2c 



(26) 



The integral in Eq. ( 26 ) can be omitted in the case if the 



z-dependence of Ap along the classical trajectory of the 
particle in the laser field is negligible. For the total phase 
of the wave function, we therefore find 

S T (x,x ) = S L (x,x ) + S P (x,x ) 

= p-{x-x ) (27) 

A(rj,x L (x,Tj,77)) 



dr/ 



'/() 



A(r7,x L (x,r7,r7)) 
t 3 ' 2c 

kh -p 



where 770 and Xq are arbitrary constants. We derive the 
propagator as follows 



G T {x,x') ps 



cd 3 q 
2e q (2^) : 



exp [ — i $>(x, x') 



(28) 

where $(x,x') = St(x,Xo) — St(x',xq). The latter can 
be represented as 



$(x, x') = S T (x, x') + A T (x, x') 



(29) 



with 



A T (x,x') 
1 rv 



A(77,X L (X,77,77)) 



A(?7,x L (x, Ty,?? 7 )) 
2c 



A(t7,x l (x,?7,77)) 



2c 



A(r?,x L (x, 77,77')) ' 



(30) 
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The two terms in the integrand of Eq. ( 30 1 deviate only in 



xl(x, 77,77) and xl(x, 77, 77'). Within the saddle-point ap- 
proximation applied later, 77 and 77' are the phase of recol- 
lision and ionization, respectively. Therefore, xl(x, 77,77) 
and xl (x, 77,77') differ only by the distance in space be- 
tween ionization and recollision which is zero. Thus, we 
can omit A^(x,x'). 



III. RELATIVISTIC PHASE-MATCHED X-RAY 
ASSISTED HHG 

In this section, the theory developed in Sec. [II] is em- 
ployed to investigate macroscopic harmonic emission in 
the relativistic regime for a HHG setup where an IR laser 
field of relativistic intensity is assisted by an x-ray field. 
In this scheme, the x-ray frequency u>x exceeds the bind- 
ing energy J PjX of the electron and thus delivers an initial 
momentum to the freed electron which can balance the 
subsequent drift motion. The laser alignment of the setup 
and an example of a classical trajectory that recollides 
are illustrated in Fig. [2] The weak counterpropagating 
IR field (brown) is important for a phase-matched macro- 
scopic response and can be ignored when discussing the 
process for a single atom. The HHG medium is a macro- 
scopic gas of multiply- charged ions. 



A. Current density 

We begin with adapting the general equation for the 
current density Eq. (151 to the present setup. Presum- 



current density of Eq. ( 15 ): 
ja(x a ,w H ,n') 

= ^Jd 4 xJ dVe iWH< ik '• 

x (j - ik H )G T (z, x')^i(x')0(x' - x a ) (31) 



"(x-x a ) 



In this expression, the x-ray field enters only into the 
potential Vi(x) = 2x • (Ex + El) in the Klein-Gordon 
formalism. Ignoring tunnel ionization by the IR laser 
field which is justified in the considered setup, we drop 
the laser field in this term and approximate V\{x) ~ 
2x • E x » x • Exoe~ i( " x *~ kx ' x) . In the last step, the ex- 
ponential function with the positive argument is dropped 
because it leads to an unphysical solution of the saddle 
point equations |30) . Due to its negligible ponderomotive 
potential, Ex can be neglected for the continuum prop- 
agation of the electron and, thus, it is neglected in the 
propagator. 



The dependence of the current density Eq. (31) on 



the position of the atom x a is given by the bound wave 
functions </>(x — x a , t). To separate out phase factors 
that highly oscillate with x a , a coordinate transforma- 
tion is applied: x = x — x a . Thereafter, time integra- 
tion is transformed to an integration over the laser phase: 
77 = cu^t — kL ■ x = co^t — kL ■ x a — kL ■ x. Finally, we 
obtain an expression for the current density that can be 
evaluated within the saddle-point approximation: 



j a (x a ,w H ,n 
x exp 



di] 



dr]' / d 3 qTO J (q,?7, 77', x a 



x exp 



-i(<Sp(q,W> x a) 



ing that for the realization of phase-matching additional 
weak fields Ap(?7,x) will be required, we employ the 
Green function Gt {x,x') of Eq. ( 28 ) instead of the Volkov 
propagator G L (x,x') in the general expression for the where Tp ig the laser pulse duration and 



H — k L 



kL)x ; 



(32) 



.c 2 (q + 



A(t;,x a ) 




A i 77 d x a )_kL 
c w L 



A(7?',x a ) k L 

( qH (e q • 

c w L 



e q w£ 



^p,x - c 2 ) 



k x |x-E xo |0 



(33) 



HHG 



FIG. 2: (color online). Geometry of the HHG process for 
a collinear alignment of the x-ray and laser field. The co- 
propagating x-ray field (orange) has a frequency above the 
ionization energy to achieve drift compensation. The weak 
IR field (brown) is employed to accomplish phase-matching. 



Further, 



Sp(-p,ri,r/,x a .)= / d?7(£ : q(77,x(x a ,77,?7))-c 2 +/p^)/wL 

Jrj' 

.(34) 

with the relativistic energy of the electron in the position 
dependent laser field given by 



Eq^x) = £ q 



'cl • q 



A(t7,x )^ A(j7,x) 
2c 



(35) 
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The wave vector of the laser is Isl = n re £ — e z and 



B. Single-atom response 



n rc [ 



(36) 



is the refractive index of the plasma, with the plasma fre- 
quency uj p — ^/4ttZ~p and the ion charge number Z. This 
way, we take into account the change of the phase velocity 
caused by the free electrons but ignore pulse deformation. 
Further, we restrict ourselves to a weakly focused laser 
field which can be treated as a plane wave and ignore the 
transversal variation of the vector potential with respect 
to the propagation direction. 

The propagation direction and the frequency of the as- 
sisting x-ray field Ex have to be chosen in such a way to 
facilitate the phase-matching and to counteract the rela- 
tivistic drift effectively Let us first consider the feasibil- 
ity of phase-matching by analyzing the harmonic emis- 
sion phase of different ions in the medium via Eq. (32). 



Generally, the vector potential consists of two terms as 
in section III CI 



A(t7,x) = A L (?y) + Ap(?7,x) 



(37) 



To understand the most favorable condition for phase- 
matching, let us assume for a moment that there are 
no additional fields for facilitating quasi-phase matching 
and ignore deformation of the laser pulse. Then, we can 
omit Ap(77,x) in Eq. (37) [a modified setup where the 
additional Ap(r/,x) field is taken into account will be 
considered in Sec. Ill C . In this case, we have no x^ de- 
pendence of m J (q, r]',r}",XA) and <Sp(q, rf , T]", x^). The 
only contribution to the variation of the harmonic emis- 
sion phase along the medium comes from the exponential 
function exp(iAkr a ) in Eq. (32) determining the coher- 
ence length l co h = 7r/|Ak|, with the phase-mismatching 
wave- vector 



Ak 



(38) 



caused by the refractive index difference between the 
fields. The phase-mismatch Eq. (38) can be split up in 
two parts. One is due to the phase-mismatch between 
the harmonics and the laser and the other one due to 
the phase-mismatch between the ionizing x rays and the 
laser. Interestingly, due to the different signs both can 
partially cancel out. This fact also holds true for this 
setup in the non-relativistic regime. Note that 
and kjj are much larger than the other two terms left 
in Eq. (38) for the case ojh ^ wx- The phase mismatch 



is smallest in the case of a collinear alignment of the x- 
ray and laser fields kn || ki, (see Fig. [2]) with propagation 
in the z direction. Accordingly, only the x component of 
the spectral current density of Eq. ( 32 ) will contribute to 



the harmonic emission in the z direction. 

The drift compensation can be achieved with an appro- 
priate choice of the x-ray frequency which is considered 
in Sec. IIIIB2I 



1. Singe-atom HHG rate 

The single-atom photon emission rate per solid angle 
of the n th harmonic can be calculated from 



dw n 1 f n ^+^ dN 
dw H 



(39) 



cR 2 /-ri^L+WL 



(2^) 2 w h Tp 



dw H |EH,o(n',w H )| 2 



using Eq. ^ for EH.o(n',WH) and the density p(x a ) 
<5(x a ) which yields 



dw n _ W^H |; 

dfl ~ (2ttc 



3 |j a (0,7io; L ,n')| 2 



(40) 



where in 

Eq. J321, Tp is replaced by 



the expression for 
Tp is replaced by ^ L 



k(0,nu L ,n') given by 
to confine the emission 



to one laser period [the rate of Eq. ( 40 1 is identical to 
ones in Refs. [HI US]]. 

A typical HHG spectrum for the considered setup cal- 
culated within the saddle-point approximation is dis- 
played in Fig. [3] (a) in blue for the set of parameters 
denoted in the caption of the figure. In this paper we 
employ a zero-range potential |50| to model the binding 
potential. Analytical expressions of the matrix elements 
which appear in Eq. (33) can be found in, e.g., |17 | 151 ) . 
The ionization potential I p K = 8 a.u. is chosen large 
enough such that tunnel ionization by the strong opti- 
cal laser field does not lead to depletion of the bound 
wave function. We compare the spectrum obtained from 
the x-ray assisted setup (blue line) with the spectrum of 
a conventional HHG setup where no x-ray field is present, 
calculated either fully relativistically [T7] (dashed black) 
or within the dipole approximation (DA) |51j (gray). The 
ionization matrix elements for the conventional HHG se- 
tups was multiplied by the factor rg^TOyi P^l t° account 
for the underestimation of the tunneling rate when em- 
ploying the zero-range potential. Here, t' is the ioniza- 
tion time and «t = •v/2I P) t, with J P)t being the ionization 
potential in the case of the conventional setup. In or- 
der to have a fair comparison of the x-ray assisted setup 
with the conventional one, we have to choose the ioniza- 
tion potential of the conventional setup (ip,t) different 
from that of the x-ray assisted one (I P:X ) such that the 
ionization rates of both setups were the same. For this 
purpose, the ionization rate for the conventional setup is 
calculated from the Perelomov, Popov, Terent'ev (PPT) 
tunneling rate |53l 154] . while the single-photon ioniza- 
tion rate from the zero-range potential is derived using 
the differential photoionization cross section in nonrel- 
ativisitic, dipole approximation and approximating the 
ionized electron wave function by a plane wave (in anal- 
ogy to Ref. [55]): 



do" 



x.lph 



V 



2ttcu>x 



I(pMo>| s 



2 4 p 3 wx^z cos 2 9 
c (n 2 +p 2 )' 



-(41) 



7 




46 48 50 52 

HHG energy u H [keV] 




46 48 50 52 

HHG energy u H [keV] 



54 



FIG. 3: (color online), 
bility for Eq — 2.5 a.u. 



(a) Single-atom emission proba- 



Sa.u., ojx 



14 a.u. and 



Ex = 0.65 a.u. (blue, dark gray) and a conventional laser 
field (i?x = 0) with Eo = 2.5 a.u. and 7 Pit = 4.8 a.u. (dashed 
black) and the same configuration in the DA (gray). For 
the second configuration, I Ptt is chosen such that the average 
tunnel-ionization rate is the same as the single-photon ioniza- 
tion rate in the case before, (b) Separate HHG yields of the 
three contributing quasi-classical trajectories for the discussed 
setup [blue line in (a)]. The dotted red contribution (long 
trajectory) is suppressed because the drift for the trajectory 
is not completely compensated (as discussed in Sec. |IIIB2[ ). 
The short trajectory contributions (solid lines in dark blue, 
dark gray and orange, gray) are nearly identical. 



where 9 is the angle between the electron momentum p 
and the polarization direction z of the x-ray field. An 
integration over all emission angles yields the total cross 
section 



&S,lph 



4k x it (2uj x - Klf/ 2 



3c 



(42) 



where we used p = \/2uj x — k%. 

The main message of Fig. |3] (a) is that the relativis- 
tic drift can be fully compensated in the x-ray assisted 
HHG setup (in the case of the chosen x-ray frequency 
wx — Ip.x — 6 a.u., the gray and blue curve are of com- 
parable order, small suppression arises from the different 
spreading behavior) . The yield of the considered setup is 
much higher than that for the conventional setup (dashed 
black), the latter being suppressed by the drift. 

In the next section, we explain how the single- x-ray- 
photon ionization provides the necessary initial momen- 



tum for the electron (opposite to the IR laser propagation 
direction) to counteract the relativistic drift in the case 
when the x-rays propagates along the strong IR laser 
field. We discuss also the optimization of the applied 
x-ray frequency for the HHG process. 



2. Drift compensation and influence of x-ray frequency 



The integration in Eq. (32 1 is carried out via the 
saddle-point integration method [HI [SS]. This means 
that instead of the integration we only need to sum the in- 
tegrand over a small number of saddle points for each en- 
ergy a>H- A saddle point (17, 77', q) determines the ioniza- 
tion and recollision times and the canonical momentum 
for the electron classical trajectory leading to the har- 
monic energy under consideration. In general, they are 
complex expressing non-classical dynamics during tun- 
neling ionization. For the parameters chosen above, 3 
quantum paths (saddle points) contribute to the spec- 
trum for each energy being equivalent to three classical 
trajectories that recollide with that energy. The separate 
contributions of each quantum path to the spectrum are 
shown in Fig. 3] (b) . The two paths marked in blue and 
orange (solid lines) have nearly the same yield whereas 
the dotted red line is suppressed by several orders of mag- 
nitude. In the following we explain the reason for the dif- 
ference and discuss the influence of the x-ray frequency 
on the dynamics. 

In order to understand the number of contributing 
trajectories in Fig. [3] (b), we calculate the saddle-point 
solutions for different x-ray frequencies wx for the har- 
monic emission at 50 keV and show the ionization phase 
saddle point in Fig. [4] (a). For small initial energies 
wx-/ p ,x, two saddle points contribute to harmonic emis- 
sion as in the usual case of HHG in a laser field only. 
Both saddle points, the long (Re 772 ~ —1.345) and short 
(Re 772 ~ —1.115) trajectory, are complex [their real part 
is shown in the graph] and their HHG amplitude is very 
tiny due to the missing drift compensation which is in- 
dicated by the complex value. When increasing wx, first 
the short trajectory and then the long trajectory split 
up into two parts. These branches are called uphill and 
downhill trajectories, respectively, because their initial 
momentum component along the laser polarization is ei- 
ther positive or negative [30] . After the splitting at about 



p,x 



4 a.u. and ujy^ — L 



p,x 



8.5 a.u., the respective 
ionization phase is purely real which indicates that the 
initial momentum is sufficient to compensate the subse- 
quent relativistic drift. The short trajectory reaches drift 
compensation earlier because it spends less time in the 
continuum and, therefore, undergoes a smaller drift that 
requires compensation. 

The dashed line in Fig. [I] (a) denotes the x-ray fre- 
quency ujx that was chosen in Fig. [3] (b). The short tra- 
jectory has two contributions (blue and orange) whereas 
the long trajectory (dotted) has only one contribution 
(red). The contribution of the long trajectory is sup- 
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FIG. 4: (color online), (a) We show the ionization phase 
saddle points of the 50 keV trajectory for different x-ray fre- 
quencies which is the same as ionization time in units of ra- 
dian. The dashed line indicates the value chosen in Fig. [3] 
The dotted and solid lines belong to the long and short tra- 
jectories, respectivly. (b) displays the initial momentum di- 
rection for different initial energies ojx — / P .x [as indicated in 
a.u. next to the respective arrow] needed for the emission 
of a 50 keV photon. The uper and lower branch correspond 
to the short up- and downhill trajectories, respectivly, where 
the color indication coincides with the one in (a). Note that 
the ionization phases 772 are different for the up- and downhill 
trajectories, x and z are the propagation and polarization 
directions of the laser, respectively, (c) Differential ionization 
rate depending on the initial energy ojx 



The considered direction is in the initial momentum direction 
determined by the saddle point equations. From (b) we see 
that the z-component is approximately p ZiC = 2.9 a.u. and the 
x-component depends on ux — Ex- 



pressed by about 3 orders of magnitude compared to 
the short contributions which is visible from Fig. [3] (b) . 
This is because the long trajectory spends more time in 
the continuum and experiences a larger relativistic drift 
which cannot be fully compensated for. In this case, the 
ionization saddle point is complex leading to a damp- 
ing in the exponential function in the respective ampli- 
tude Eq. (32). By increasing u>x above u>x — Ip ~ 8.5 a.u., 



the drift compensation could also be achieved for the long 
trajectory and the dashed red contribution in the spec- 
trum could be enhanced leading to a larger single-atom 
yield. However, only one of the trajectories can be phase- 
matched in many cases and the enhancement of the other 
trajectories would not be useful. 

In Fig. [4] (b), the initial momentum vectors of the ion- 
ized electron p(r]',q) = q + - ^(e q + - c 2 ), 
which correspond to solutions of the saddle-point equa- 
tions, are displayed for different ox- When ox — I P ,* ~ 
4.2 a.u., the momentum required for drift compensation 
of the short trajectories is just reached [see Fig. [4] (a)]. In 
this case, the initial momentum is directed mainly along 
the z-direction [arrows marked with 4.2 in Fig. [4] (b)]. 
When u>x is increased, only the p x component changes; 
the p z component approximately remains constant be- 
cause it is determined by the drift compensation condi- 
tion. The electrons with an appropriate initial momen- 
tum vector can be provided by the x-ray single-photon 
photoionization because the latter happens with a large 
angular distribution with a maximum around x-ray po- 
larization direction as can be seen from the differential 
ionization cross section of Eq. (41). Because the HHG 
amplitude for each trajectory contains the differential 
ionization cross section Eq. (41), the efficiency in each 



case depends on the scalar product between required ion- 
ization direction p and x-ray field polarization direction. 
This results in some freedom in choosing the direction 
of Ex- Only if p and Ex were close to perpendicular 
[9 f=a 7r/2 in Eq. pi])], the differential ionization prob- 
ability would be close to zero. For realization of phase- 



matching, as it is shown above in Sec. Ill A the collinear 
propagation of the laser and x-ray field is advantageous. 
For this case of a collinear alignment, we show the differ- 
ential ionization probability drr x ,iph/df2 from Eq. (41 ) for 
different x-ray frequencies wx in Fig. ||(c). The emission 
angle of interest is estimated by the initial momentum 
p = y/2{uix — I P x ) and its z-component p ZjC = 2.9 a.u. 
taken from Fig. [4] (b) via sinf9 = p ZyC /p. For initial en- 
ergies just above 4 [e.g., 4.2 corresponding to the nearly 
horizontal vectors in Fig. [4] (b)], dcr x lp h/df2 is vanish- 
ing because the momentum direction and the direction 
of Ex are perpendicular. The angle of emission 8 will 
increase with rising u>x, increasing the ionization prob- 
ability. On the other side, large values for w x decrease 
the overall ionization probability due to the denominator 
of Eq. (41 ) . These two competing tendencies creates the 
maximum in the ionization probability in Fig. [4] (c) . We 



see that the chosen value wx — I] 
close to the optimal conditions. 



p,x 



6 a.u. (wx = 14) is 
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The former trajectory-based discussion can also be 
seen from a wave packet perspective. The single-photon 
ionization mechanism with a large initial kinetic energy 
obeys a dipole angular distribution of the ejected wave 
packet, i.e. the wave packet has an increased spreading 
velocity compared to tunnel ionization. The increased 
spatial dimension of the recolliding wave packet is ex- 
ploited to overcome the drift. 



C. Macroscopic HHG emission 

After discussing the single-atom yield of the x-ray as- 
sisted setup, we continue to elaborate on the macroscopic 
aspect of the emission from a gas target. 
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FIG. 5: (color online), (a) Real part of spectral component of 
the locally emitted HHG field at 48.6 keV at different positions 
along the propagation direction. The blue dashed line is for 
HHG without the quasi-phase-matching scheme. The red line 
is for the case of adding the weak counterpropagating field 
to achieve QPM. (b) The macroscopically emitted spectral 
photon number via Eq. |9| is displayed for the QPM scenario. 



We inspect the emission from a Be 3+ gas of homoge- 
neous density p = 5 x 10 16 /cm 3 with the same parame- 
ters as in Fig. [3] The plasma refractive index at the laser 
frequency is — 5 x 1CP 5 . The phase mismatching 



wave- vector of Eq. ( 38 ) at the harmonic emission energy 
of 50keV is then Ak — 6 x 10~ 4 a.u. and the coherence 
length l co ^ = -k / Ak = 0.25 /im. In order to increase 
the coherence length, a quasi-phase matching (QPM) 
scheme can be employed [45H47] , We propose to use a 
weak counterpropagating IR field with the parameters 
E 2 = 5 x 10~ 5 a.u. and lo 2 = 0.0418 a.u. to achieve quasi- 
phase matching (QPM). The additional field is denoted 



by a brown line in Fig. [2] It is included into our mathe- 
matical formalism by Ap(n, x) in Eq. (37). In that way a 
dependence on x a is introduced into 5p(q, n, 7/, x a ) and 



the saddle points for the integration in Eq. ( 32 1 depend 
on the position within the medium. Thereby, x(x a , 77, 77) 
in Eq. ( 34 ) contains the variation of the weak field seen by 
the electron along the z-direction. For the chosen set of 
parameters, the approximation of x(x a ,7y, 77) by x a does 
not lead to a significant change of the final results and, 
thus, can be done to save computation time. The impact 
of the additional field can be observed in Fig. | 5| (a ) . The 
real part of the spectral current density Eq. ( 32 1 at the 
respective position is shown. The emitted total field is 
given by a spatial integral over all contributions of the 
current density [see Eq. Q] . Without QPM (blue dashed 
line) , the single-atom contributions oscillate on the scale 
of the coherence length estimated previously. An inte- 
gration over all contributions results in extensive can- 
cellation. However, when applying the additional field, 
the symmetry between the positive and negative contri- 
butions is broken (see the red line) and both parts only 
partially cancel thus achieving quasi-phase-matching and 
a nonzero value of the integral. The parameters of the 
additional field were chosen to optimize the photon en- 
ergy at 48.6 keV. A medium length of 100 /mi was chosen 
whereas the diameter is 500 /im. The assumed laser and 
x-ray pulse duration is 10 cycles. The length is limited 
due to the assumed bandwidth of the weak QPM field 
A102 ~ 0.1%. The spectrum is shown in Fig. |] (b) and 
an integral over the spectrum yields the final result of 
5 x 10~ 7 emitted photons per shot. The number is of 
similar order of magnitude as in the other relativistic 
HHG setup based on the counterpropagating attosecond 
pulse trains for driving harmonics |25| . 



D. Efficiency analysis 

We continue with a discussion about the small HHG 
yield in the relativistic regime and identify several rea- 
sons for it that are either general to the relativistic regime 
or specific to this setup. First, we specify an estimate ex- 
pression for the emitted photon number: 



N = 



x An x AQ x At x V 2 p 2 



(43) 



It allows in a simple way to estimate the HHG yield by 
an order of magnitude and to single out the different is- 
sues influencing the HHG yield. In Eq. (43), dw n /dfl is 



the single-atom emission rate, An is the number of har- 
monics within the phase-matched frequency bandwidth, 
At the interaction time that is approximately the delay 
between both pulses, Ail the solid angle of emitted har- 
monics, V the volume of coherently emitting atoms [per- 
fect phase-matching is assumed in this volume], p the 
atomic density. 

First, we demonstrate the usefulness of the expres- 
sion by estimating the photon number for the proposed 



10 



setup and show that the result of the former exact cal- 
culation can be reproduced. We estimate the terms in 
Eq. (43 1 as follows: the single atom emission rate is 
dWn/dTl « 10~ 21 [see Fig. [3] (b)]; the phase- matched 
frequency bandwidth can be deduced from Fig. [5] (b): 
Awh ~ 0.2 keV which gives An ~ 10 2 ; the solid an- 
gle of phase-matched emission in the far field is deter- 
mined by the interference pattern of a circular aperture 
AQ ~ tt{2ttc/Duj b ) 2 w 10~ 14 , where w H = 50keV and 
medium radius r a = 5 x 10 6 a.u. are assumed; the in- 
teraction time At ps 10 3 a.u. is the laser pulse duration; 
the volume is cylindrical V = 7rr 2 Az a = 10 20 ; the plasma 
density restricted by dispersion is p = 5 x 10 /cm 3 . Tak- 
ing all pieces together, the emitted photon number under 
phase-matched conditions is 



a r50 kcV 
iv rel 



dw 7 , 



dfl 

= 10~ 21 x 10 
= 10~ 6 



x An x AS] x At x V 2 p 2 

x,s 



(44) 

x 10~ 14 x 10 3 x 10 40 x 10~ 16 



in agreement with the previous accurate calculation. The 
subindex x stands for the x-ray assisted setup whereas s 
denotes that the short trajectory contribution is taken 
into account only. 

To explain the low yield in the relativistic regime we 
perform the same kind of estimation for a state-of-the- 
art HHG experiment [57] in the non-relativistic regime 
where 80 as pulses were generated from harmonics below 
100 eV. We can estimate the emitted photon number in 
this case: 



N, 



50 cV 



dw r . 



1 — lcl 



dSl 



2 „2 



= 10" 

= 10 9 



-16 



x An x AO x At x V z p 
x 10 1 x 10~ 7 x 10 2 x 10 41 x 10" 



(45) 



The single-atom contribution dw n /d£l\ t , s was calculated 
from |51) where we additionally inserted a correction fac- 
tor accounting for the underestimation of the tunneling 
rate when using the zero-range potential as described 
in [55] . By comparing Eq. ( 44 1 and Eq. ( 45 ) , one observes 



a dramatic suppression of 15 orders of magnitude when 
rising the HHG energy by about 3 orders of magnitude. 
It arises mainly due to the single-atom yield dw n /dQ, 
the phase-matched emission angle Af2 and the gas den- 
sity. The single-atom contribution will be investigated 
separately below. The estimated solid angle emission 
angle decreases quadratically with the harmonic energy. 
This is because a smaller harmonic wavelength leads to a 
smaller angle of the first interference minimum. The gas 
density depends on the phase-matching conditions which 
are much more difficult to fulfill in the relativistic regime 
and thus the gas density is lower in this case. 

In the following, we inspect the ratio between the 
single- atom yields dw n /dCl of Eq. (44 1 and Eq. (45) 



closer. In each case, we concentrate on a single (short) 
trajectory at 50 eV and 50keV, respectively. Then each 



single-atom rate can be estimated [T5J [TU] , 



dw r . 



dn 



1 



(2ttc) 3 w 



(46) 



P(«,q)> (p(7/,q)|£ X z|0) 



for the considered relativistic setup [see also Eq. (40)], 
and by [5T] 



dw n 



dfl 



(2ttc) 



3^H 



1 



{-2-Kif 



D(p,r),r]') 



A(n)/c) 



(47) 



2\/2k 3 



\E(V')\ 



^-(p ;c +A(n')/c|F|0)e- i(S(q ' f '' f '' )+ ^ i 



for a conventional nonrelativistic setup where we inserted 
the tunneling correction factor of Ref. [52] and where 

D(p,r),r]') 



det 9 ^ p ' v ' v ,? is the functional determinant 



of the respective action. Both, Eq. (46) and Eq. (47) 



are evaluated at the saddle point belonging to the short 
trajectory of the respective energy. All factors in Eq. (46 ) 



and Eq. ( 47 1 were ordered in the same way and a distinct 



physical meaning can be assigned to them [251 EE] : 



dw n 
~dO~ 



P/l 



dti 



= ( Pf )| 2 ^%^. (48) 
dpi 3 



The factor accounts for the phase space and converts 
the matrix element into the probability, d 3 u>i(ii, pi)/dpi 3 
is the differential ionization rate with the ionization time 
ti in momentum direction pj, a rcc (Pf) is the recombina- 
tion amplitude and the last factor accounts for the dy- 
namical properties of the wave, pf is the final momentum 
at recollision, duju/dti is the electron wave packet chirp- 
ing factor discussed in |25j . 

We compare all factors in Eq. (48) between the two 



cases, to identify the reasons for the five orders of mag- 
nitude suppression of the single-atom yield in the rela- 
tivistic regime. Since the harmonic energy increases by a 
factor of 10 3 by going to relativistic case, the factor us^ 
yields an increase of 6 orders of magnitude. The differ- 
ential ionization rate of the particular trajectory is re- 
duced by a factor of 10~ 2 . Three properties contribute 
to estimate this ratio: The electron angular distribution 
of ionization is much broader for the one-photon ioniza- 
~ -^p.x) ~ 6 than for tunneling ionization 
0.23 [19] yielding a factor of 4 x 10~ 2 . 



On the other hand, the total (constant) ionization rate 
of the relativistic example [see Eq. ( |42[ )] is by a factor 
2 higher than the instantaneous rate of non-relativistic 
example 53, 54 . Third, in the relativistic setup, the rel- 
evant electron trajectory starts with a certain angle 9 off 
the x-ray field direction resulting in a factor cos 2 9 ss 0.3 
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in Eq. (42). Together, we find the ratio of the differen- 
tial rate 0.3 x 2 x 4 x 1CT 2 « 1CT 2 . Note that the total 
time-average ionization rate is a factor of 10 lower for the 
non-relativistic example than for the relativistic one [see 
Eq. (42)]. The recombination amplitude is reduced by a 
factor of 10~ 5 as discussed in [25]. The factors pf and 
duiY{/dti contained in the functional determinant reduce 
the relativistic yield by 10 -4 . Collecting all factors, we 
find the suppression of 10 6 x 10~ 2 x 10" 5 x 10~ 4 = 10~ 5 
of the single-atom yield according to ratio between the 
respective terms in Eq. ( 44 ) and Eq. ( 45 ) . 



IV. CONCLUSION 

Extending table-top HHG to the hard x-ray domain is 
an exciting prospect, especially because many research 
labs already use HHG as a XUV sources and other ap- 
proaches to generate hard x rays, like free electron lasers, 
require large scale facilities. 

The present study discussed several difficulties that 
need to be overcome in order to realize the idea. The 
relativistic drift has been extensively discussed in the lit- 
erature. Each proposed geometry has its own advantages 
and disadvantages regarding phase-matching. Generally, 
increasing the harmonic energy renders phase-matching 
more difficult for many reasons: the emission phase of 
the harmonics depends approximately linearly on the in- 
tensity (</> ~ U p t). Small intensity variation, e.g., in a 
Gaussian focus, immediately results in phase difference 
much larger than it. On the other hand, differences in 
the phase velocities between the harmonics and the laser 
wave lead to a slip in space between both waves. This 
results in phase-mismatch as soon as the slip is compara- 
ble to the harmonic wavelength which happens earlier for 
shorter harmonic wavelengths. Additionally, relativistic 
HHG is always accompanied by a large ionization lead- 



ing to an enormous plasma dispersion. For these reasons, 
in the best case, we obtain realizable medium lengths of 
only a few tens of urn reducing the expectable macro- 
scopic yield. 

Apart from the relativistic drift and phase-matching, 
we identified further issues decreasing the harmonic emis- 
sion in relativistic HHG connected with the single-atom 
yield. First, recombination of the recolliding electron be- 
comes less likely for high momenta: scattering is favored 
instead. Secondly, the electronic wave function is spread 
over a larger energy bandwidth. If phase-matching can- 
not be achieved for the whole bandwidth, however, a 
large part of the harmonic radiation is lost. This was ex- 
pressed by the chirping factor. Third, the solid emission 
angle decreases quadratically with the harmonic energy 
increase. 

Regarding the harmonic yield, the present setup for 
relativistic HHG as well as the one in [25] yield a small 
photon number for emitted harmonics that are both of 
similar order. On the bottom line, we think that the 
setup considered in this paper is more promising than 
that of [25 because the required laser intensities are lower 
and the phase-matching scheme is more practical. 

One important conclusion of the paper is that phase- 
matching favors the collinear alignment of the laser and 
x-ray beams for the x-ray assisted relativistic HHG setup. 
This co-propagation is sufficient to induce drift compen- 
sation and no perpendicular alignment of both beams is 
required. Note that the collinear geometry has already 
been used in various experiments [26, 2 71 129] . 
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